Electric field switching of antiferromagnetic domains in YMn 2 Os: a probe of the 
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We employ neutron spherical polarimetry to determine the nature and population of the coexisting 
antiferromagnetic domains in multiferroic YM^Os. By applying an electric field, we prove that 
reversing the electrical polarization results in the population inversion of two types of in-plane 
domains, related to each other by inversion. Our results are completely consistent with the exchange 
striction mechanism of ferroelectricity, and support a unified model where cycloidal ordering is 
induced by coupling to the main magnetic order parameter. 

PACS numbers: 25.40.Dn, 75.25.+Z, 77.80.-e 



For a long time, ferroelectricity and magnetism have 
been thought of as essentially independent. Recently 
however, a class of "novel" multiferroics, such as 
TbMn0 3 Q and REMn 2 5 (RE=Y, rare earth, Bi) 
has been discovered, in which antiferromagnetic or- 
dering itself is though to induce ferroelectricity. Unlike 
'conventional' multiferroics, magneto-electric coupling in 
these materials is very strong, leading to spectacular ef- 
fects, such as rotation or reversal of the electrical po- 
larization upon application of a magnetic field The 
mechanism (or mechanisms) leading to ferroelectricity 
upon magnetic ordering are not known with certainty [J] . 
Most theories 0, Q have been guided by the observation 
that the majority of these materials possess a cycloidal 
structure, in which the spin direction is rotated contin- 
uously, and have postulated that ferroelectricity arises 
from displacements, due to spin-orbit coupling, of ionic 
or electronic density between two metal ions whose spins 
are non-collinear. In its high-temperature commensu- 
rate phase. The YM^Os magnetic structure is almost 
coUinear [3, @] - its spins are only slightly tilted ( 15°) 
away from collinearity - but it has been shown [8[ that 
this small tilt, including a weak cycloidal component in 
the crystallographic be plane, could in principle induce 
ferroelectric displacements through spin-orbit coupling. 
An alternative mechanism, not requiring non-collinearity, 
has been proposed Q for YMn 2 05, in which small dis- 
placements of oxygen or metal ions optimize magnetic su- 
perexchange energy (exchange striction). Regardless of 
the specific mechanism, a consistent prediction for these 
'novel' multiferroics is that different directions of the elec- 
trical polarization should be associated with distinct an- 
tiferromagnetic domains. In general, one would expect 



these domains to be related by the symmetry operators 
that are lost at the ferroelectric transition, namely the 
center of inversion. However, the case of REMn.205 is 
not trivial in this respect: if only one if the two com- 
ponents (cycloidal or planar) is coupled to the electrical 
polarization, it is conceivable that only that component 
will respond to an applied electric field. This possibility is 
particularly appealing in the case of the spin-orbit mecha- 
nism, since in this case the polarization could be reversed 
by switching the small c-axis component S z , while leav- 
ing the large in-plane component S xy , and consequently 
the main antiferromagnetic domain structure, unaltered. 
Whether this will happen depends on the strength of the 
magnetic coupling between the two components, which, 
however, can be very weak, since the only coupling terms 
allowed are antisymmetric Observation of these do- 
mains and of the domain switching by an external elec- 
tric field provides therefore a powerful insight in the cou- 
pling of the two substructures to the polarization and to 
each other. Optical second-harmonic generation (SHG) 
ficl 11 1 and scattering of polarized neutrons [T3j have 
been previously employed to observe domain formation 
and switching in multiferroics. However, SHG is intrin- 
sically limited to magnetic ordering at the gamma point, 
and scattering of polarised neutrons (without analysis) 
affords limited information about the structure of the do- 
mains, which must be known in advance. In this Letter, 
we describe a direct observation of the YM112O5 antifer- 
romagnetic domains using neutron spherical polarime- 
try, which is uniquely sensitive to both domain structure 
and domain population. By cooling the crystal in posi- 
tive or negative external electric fields, we determine the 
in-plane magnetic structures of the two domain types 
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uniquely, and conclude that switching from one to the 
other involves reversal of the in-plane components S xy of 
the spins, not just a change in the small S z component. 
The results of this experiment are most easily under- 
stood in the framework of the exchange striction model, 
where S xy is directly coupled to the polarization, and also 
suggest that S z is induced by the direct Dzyaloshinskii- 
Moriya (DM) interaction. 

The YMn 2 5 magnetic structure for 24 K < T < 38 K 
has been determined from neutron powder [7|] and single 
crystal [H,[l3| data, and is shown in Fig Q] for the possible 
polar domain configurations. It consists of staggered an- 
tiferromagnetic zig-zag chains running along the crystal- 
lographic a axis (horizontal in the figure, panels I and II), 
with spins parallel within each chain and tilted by about 
15° with respect to the a axis. Configurations I and II are 
related by inversion, although they can also be obtained 
by reversing the direction of the spins in half of the chains 
[131 ] . The sign of the out-of-phase c-axis component (con- 
figurations III and IV also related by inversion) deter- 
mines the rotation direction of the cycloidal modulation 
in the frc-plane. Unique domains are obtained by combin- 
ing in-plane and cycloidal components: the combination 
I+III is related by overall inversion symmetry to II+IV 
and so is I+IV with II+III. Domains differing solely by 
the sign of the c-axis component, like I+III and I+IV 
are not related by symmetry and are in principle distin- 
guishable by diffraction. Indeed, neutron diffraction data 
clearly favor the I+III and II+IV configurations. How- 
ever, an admixture with the other domains would result 
in a reduction of the refined S z component, and can not 
be completely ruled out. As already mentioned, the en- 
ergetic degeneracy between these domains is lifted only 
by weak antisymmetric exchange |{|. 

The single crystals of YM112O5 was grown using 
B 2 3 /PbO/Pb flux in a Pt crucible. The flux was held at 
1,280 °C for 15 hours and slowly cooled down to 950 °C 
at a rate of 1 °C per hour. Crystals grew in the form of 
cubes. The sample dimension was 4x0.8x4 mm 3 along the 
crystallographic a, b and c directions. Thin silver-paint 
electrodes were attached to the crystal faces perpendic- 
ular to the b-axis, with gold wires providing electrical 
connections to the voltage supply. The sample was at- 
tached to an Al support with " GE" varnish, and mounted 
with the b-axis perpendicular to the scattering plane. 
Neutron spherical polarimetry data were collected as a 
function of temperature using the diffractometer D3 at 
the Institut Laue-Langevin (Grenoble, France) equipped 
with CRYOPAD-II. The crystal was mounted with the 
b-axis (i.e., the direction of the electrical polarization) 
perpendicular to the diffraction plane (vertical), so that 
only ZiOZ-type Bragg peaks were measured. An external 
electric field of up to 2.2 kV/cm and of either polarity was 
applied to the crystal along the b axis. Neutron spheri- 
cal polarimetry is described in detail elsewhere b4| and 
also summarized in the EPAPS supplementary material. 



Briefly, we can align the spins of the incident neutrons in 
any chosen direction and determine both the magnitude 
and direction of the polarisation of the scattered beam. 
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FIG. 1: (color online) left: schematic representation of the 
neutron spherical polarimetry experiment for the two do- 
mains, here shown in the idealized case of an unpolarized 
incident beam. The real (or imaginary) component of the 
magnetic structure factor projection M^ fci rotate clockwise 
(counter clockwise) for Domain I (Domain II), creating a 
spin polarization of opposite signs for the scattered neutron. 
The direction of the applied electric field is also indicated, 
right: Magnetic structures of YM^Os for the different do- 
main configurations, projected on the ab plane (I and II) and 
the be plane (III and IV). Small arrows represent magnetic 
moments. The observed domain switching mechanism is rep- 
resented by the inversion (change from the red to the blue) in 
the central chain (between configurations I and II in the ab 
plane). 



The polarization of the scattered beam for each Bragg 
peak hkl can be calculated for any incident polarization 
from the nuclear structure factor Nhu and the magnetic 

" [3. Like the 
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structure factor M.hki of the crystal 
familiar X-ray structure factor, N^u is a complex num- 
ber, whereas M^; is a complex vector. Neutron diffrac- 
tion and polarimetry are only sensitive to the projec- 
tion of M/,fci perpendicular to the scattering vector Q /lfeJ , 
here indicated as M^ fei . The equations for neutron po- 
larimetry are greatly simplified when nuclear and mag- 
netic scattering do not interfere and when the incident 
beam is fully polarized, as in the present case 17| . These 
equations are best expressed in the so-called Blume refer- 
ence frame [15| , with the X-axis parallel to the scattering 



vector (i.e., bisecting the incident and scattered beam), 
the Z-axis perpendicular to the scattering plane (vertical 
in our case) and the Y-axis completing the right-handed 
set (see Fig[T|). Noting that M ftj y lies in the YZ plane, 
without loss of generality we can write: 
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Where y and z are unit vectors in the Y and Z direc- 
tions. If the incident beam is fully polarized, the equa- 
tions are further simplified, and can be expressed using 
the matrix elements Py (i,j=K, Y , Z), which represent 
the polarization measured in the direction j if the inci- 
dent beam is polarized in the direction i 
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These matrix elements result from two terms: the ro- 
tated polarization, which is proportional to the incident 
polarization, and the created polarization (elements P yx 
and P zx ), which is present even if the incident beam is 
unpolarized, and is always parallel to the scattering vec- 
tor (i.e., to X). The best way to visualize Eq. Q] is the 
following: for all but the simplest magnetic structures, 
the phase factor ip in [T] varies from one chemical unit cell 
to the next along the so-called "propagation" direction, 
so that the real and imaginary components of M^" fc; de- 
scribe an ellipse in the YZ plane (FigfTJ) . The parameters 
a and <fi define the cardinal equation of the ellipse, while 
the sign of </> establishes whether the rotation around the 
ellipse is clockwise or counter clockwise. The unique sen- 
sitivity of spherical polarimetry to the domain structure 
stems from the fact that for the same reflection hkl the 
rotation is opposite for inversion domains, and conse- 
quently the signs of the created polarization P yx and P zx 
are reversed (Fig. [I]). Crucially, the magnetic structure 
itself need not be chiral or even rotating, but cannot be 
perfectly collinear, because in this case the created po- 
larization will vanish. 

In the geometry we employed, the polarization matrix 
elements are essentially insensitive to the small c-axis 
spin component [H, Ha ], so we only probe domain con- 
figurations I and II, which, as already mentioned, are 
related by inversion. The parameters and for this struc- 
ture are readily calculated for each Bragg reflection hkl. 
In particular, for Domain I </> w +90° for h = 2n + -|, 
<f> k —90° for h = (2n + 1) + |, while the signs are re- 
versed for Domain II. The parameter a depends on hkl, 
but does not change between the two domains for a given 
Bragg peak. 

The results of our spherical polarimetry measurements 
at 25 K are displayed in Fig. [2] The main panel shows the 
significant matrix elements for the reflection — \ — | , 
while the inset summarizes all results for several Bragg 
peaks. Data on the left and right panels were collected 




after cooling the crystal through the magnetic transition 
(Tat =40 K) under a positive or negative electric field, 
respectively (E= ± 2.2 kV/cm). The calculated values, 
based on the previously described magnetic structure and 
either 100 % population of Domain I or Domain II (no 
adjustable parameters), are in excellent agreement with 
the observations. In particular, the sign of the created 
polarisation (P yx = P zx ) is reversed between the two 
electric field polarities, whereas the magnitudes and the 
signs of the other elements are unchanged, as expected 
for a complete population reversal. 
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FIG. 2: (main: observed (solid/color and error bars) and 
calculated (rectangles) neutron spherical polarimetry matrix 
elements at 25 K for the — | — | Bragg peak of an YM^Os 
crystal cooled in a positive (left, red color) and negative (right, 
blue color) electric field of ± 2.2 kV/cm. Note how the signs 
of the created neutron polarization elements P yx and P zx are 
reversed between the two field orientation, while all the abso- 
lute values and the signs of the other elements stay the same 
- a clear indication of domain population reversal, insets 
observed and calculated matrix elements for several Bragg 
peaks in the two field orientations. The model has no ad- 
justable parameter, and assumes 100 % population of one of 
the two domains. (color online). 

At 25 K, the domains are strongly pinned, and cannot 
be switched by reversing the electric field at constant 
temperature. Therefore, we warmed the crystal to 35 K, 
which is sufficiently close to the transition to observe field 
switching. The domain population was still biased by the 
previous field cooling, so we could only observe a partial 
hysteresis loop (Fig. [3|). Here, the sign of the created 
neutron polarization is reversed, but the magnitudes are 
not the same for +2.2 kV/cm and -2.2 kV/cm, indicative 
of domain ratios 80 %-20 % and 43 %-57 %, respectively. 
In a separate experiment, the electrical polarization of 
a crystal of the same batch was measured under similar 
conditions. Both the hysteresis bias and the switching 
ratios were found to be in very good agreement with our 
neutron data (inset to Fig [3]) 

The simplest version of the spin-orbit model, directly 
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FIG. 3: main: partial hysteresis loop measured on the cre- 
ated neutron polarization element P yx for the — i — ? Bragg 
peak of an YM^Os crystal, warmed to 35 K after previous 
cooling to 25 under a negative -2.2 kV/cm electric field, in- 
set: integrated pyroelectric currents measured on a 0.5 mm 
thick YMn2C>5 crystal of the same batch on cooling to 25 K 
in a negative -2.0 kV/cm electric field (bottom/blue curve), 
followed by warming to 35 K and switching to a positive +2.0 
kV/cm electric field (top/red curve). The data are normal- 
ized to the fully saturated value at 25 K. Both hysteresis bias 
and the switching ratios are in very good agreement with the 
neutron data, (color online) 

derived from TbMn03 and related materials, combined 
with the observation that the magnetic coupling between 
cycloids and planar structure can be weak, naturally 
leads to the prediction that only S z should reverse with 
the electric field — a prediction that is strongly contra- 
dicted by our data. Conversely, the reversal of the in- 
plane spin components is a strong prediction of the ex- 
change striction mechanism (18(, and is here completely 
verified. We remark, however, that our data are com- 
pletely consistent with a electric-field-driven transition 
between inversion-related domains in which both the S xy 
and S z components switch in alternate chains, fn this 
scenario, which we believe to be the most plausible, the 
in-plane and cycloidal components would always have 
the same polarity for a given direction of the polariza- 
tion, and could in principle both contribute to ferro- 
electricity. The relative importance of the two effects 
can however be gauged based on the same neutron data 
for YMJI2O5, the cross product Si x Sj (related 
to the spin-orbit mechanism) is about 40 times smaller 
than the dot product Si • Sj and a factor of 100 smaller 
than Si x Sj for multiferroic TbMnOa which displays 



a spontaneous polarization similar to that observed for 
YM112O5. Therefore, switching of the S z component, if 
observed, would indicate that antisymmetric coupling be- 
tween cycloids and planar structure is important for the 
overall stability of the magnetic structure, strongly sug- 
gesting that the spirals are in fact stabilized by direct 
Dzyaloshinskii-Moriya (DM) rather than next-nearest- 
neighbor interaction. A unified picture would therefore 
emerge, in which the frustrated geometry of the Mn ions 
is primarily responsible for breaking inversion symmetry 
upon magnetic ordering, while both cycloidal ordering 
and ferroelectricity are induced by coupling to the main 
polar order parameter. 
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